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Abstract 

Effects of neutrino free streaming is evaluated on the primordial spectrum of gravitational radi- 
ation taking both neutrino chemical potential and masses into account. The former or the lepton 
asymmetry induces two competitive effects, namely, to increase anisotropic pressure, which damps 
the gravitational wave more, and to delay the matter-radiation equality time, which reduces the 
damping. The latter effect is more prominent and a large lepton asymmetry would reduce the 
damping. We may thereby be able to measure the magnitude of lepton asymmetry from the 
primordial gravitational wave spectrum. 
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I. INTRODUCTION 



Cosmological baryon asymmetry and lepton asymmetry are among the most fundamen- 
tal quantities of our Universe. The magnitude of baryon asymmetry is determined by two 
independent methods. Baryon asymmetry can be estimated by comparing the prediction of 
abundances of the hght elements with observations of those quantities in big bang nucle- 
osynthesis (BBN) [1]. The observations of small scale anisotropies of the cosmic microwave 
background (CMB) also determine the baryon asymmetry. Recent results of the Wilkinson 
Microwave Anisotropy Probe (WMAP) coincide with the results of BBN within two sigma 
level [21, 0]. Thus, cosmological baryon asymmetry has been determined very precisely. 

On the other hand, lepton asymmetry has still large uncertainties both theoretically and 
observationally. At very high temperature, sphaleron processes are active, to convert lepton 
asymmetry to baryon asymmetry and vice versa Then, it is naively expected that the 
magnitudes of baryon and lepton asymmetry are of the same order through these effects. 
But, it does not always hold true and many scenarios which allow large lepton asymmetry 
have been proposed 0, 0, Isl, |9|. Then, the determination of the magnitude of lepton 
asymmetry is crucial to determine cosmic history of our universe. Contrary to the case of 
baryon asymmetry, however, cosmological lepton asymmetry cannot be observed directly 
and only upper bounds on lepton asymmetry has been obtained by indirect methods such 
as BBN, CMB, large scale structure, and so on. The most stringent constraint comes from 
BBN and the degeneracy parameter ^ is constrained as 

|e.J<0.2 or |e.,,..|<2.6 (1) 



for the case of no mixing 10| ^ and 

lel < 0.07 (2) 



for each family in the case of strong mixing Here the degeneracy parameter ^ is the 
ratio of the chemical potential fi to the (effective) temperature, 

^ = jr- (3) 

In this paper, we argue that lepton asymmetry may be directly measurable if one can 
observationally determine the primordial gravitational wave spectrum precisely. During 
inflation, stochastic gravitational waves are generated and expanded to cosmological scales 
[l3| . Then, such gravitational waves memorize all cosmic history during and after inflation, 
and provide us with fruitful information. For example, the amplitude of the spectrum gives 
us the energy scale of inflation directly. The spectral shape reflects the history of cosmic 
expansion: the spectrum is roughly proportional to /~^, /°, / for the modes which reenter 
the horizon in the matter dominated phase, the radiation dominated phase, and the kinetic- 
energy dominated phase. This feature may be used to probe the possible change of the 
equation of state in the early universe Recently, Watanabe and Komatsu took into 

account the change of the effective number of degrees of freedom of all standard model 
particles of particle physics and showed that such changes leave characteristic features in 



^ In case that the solution of the solar neutrino problem is other than the large mixing angle MSW solution 
and that the angle ^13 is sufficiently small but non-zero, flavor mixing can be suppresse d J 1.1 . Furthermore, 
the introduction of hypothetical coupling of neutrino can also suppress flavor mixing [I4I . 



2 



the spectrum [l5|.^ Furthermore, Weinberg pointed out that the square amphtude of the 
gravitational wave is reduced by 35.6% if neutrino free streaming effects are taken into 
account |17l] . 

In this paper, we evaluate neutrino free streaming effects on the primordial gravitational 
wave spectrum by taking into account lepton asymmetry (chemical potential of neutrinos) 
and neutrino masses. Lepton asymmetry induces two opposite effects. One is to increase 
anisotropic pressure, which increases the damping of the spectrum by free streaming effects 
of neutrinos. The other is to delay the matter-radiation equality, which decreases it. Then, 
we show that net effect of lepton asymmetry is to decrease the damping of the amplitude of 
the primordial gravitational wave by free streaming effects of neutrinos. 

In the next section, we briefly review the basics of lepton asymmetry and derive 
anisotropic pressure of neutrino in the presence of lepton asymmetry. In Sec. Ill, we 
investigate the evolution of gravitational wave and discuss the damping of the amplitude of 
the primordial gravitational wave by free streaming effects of neutrinos. Finally we give our 
discussion and summary. 



II. LEPTON ASYMMETRIC COSMOLOGY 



In this section we quickly review the basics of massive degenerate neutrinos. Obser- 
vational implications of degenerate massive neutrinos for density perturbations have been 

3- 



discussed in the literature [18|, [19 
gravitational wave background. 



Here we apply them to the evolution of primordial 



A. background 

The ei ierg y density of one species of massive degenerate neutrinos and anti-neutrinos is 
given by [l8|, l21[ 



= {ksT^f j 0^s'im + fM) ■ (4) 

Here T^, and Tj^^ are the temperature of neutrinos and that of today, respectively, a is a 
cosmic scale factor, q is comoving momentum in units of ksT^^, m is normalized neutrino 
mass defined by 

and e = a/?^ + a^m^ is the comoving energy density. Note that comoving momentum q 
is related to the proper momentum p hy p = q/a, which redshifts as a~^. Therefore, q is 
constant relative to the expansion of the universe. When the universe was dense and hot 
enough, neutrinos were kept in thermal equilibrium with the rest of the plasma through weak 
interactions so that the distribution function had been relaxed to the form of Fermi-Dirac 



The effect of quark gluon plasma phase transition is discussed in 
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distribution. Even after neutrino dccoTipling at ~MeV, owing to Liouville's theorem, the 
distribution function is still given by the distribution: 



1 



1 + ei-^ ' 



Here ^ is the degeneracy parameter. 
The pressure is given by 



HP.+p.)^{kBT.rl (/.(?)+/.(?)) . 



1 



1 + ei+^ 



(6) 



(7) 



Here and hereafter, we assume that fi (and C,) is spatially homogeneous for simplicity. In 
our calculation we also assume, for simplicity, that all the three generations have the same 
properties, that is, we set m,,^ = m^^ = m^^ = and = = = ^. 

It will be useful to note the expressions for the energy density and the pressure of massive 
degenerate neutrinos, in the relativistic and non-relativistic limits for numerical stability and 
efficiency. In the relativistic limit, the energy density can be written as 



Pu + Pp 



27r2 

{Puo + PPo 



q'dqif, + /p) + / q'dqiU + h 




Here we have used the well-known expression for the energy density of massless non- 
degenerate neutrinos 



Pvo + pDo = {ksTj^y 



(Pq 2 
1<1- 



(27r)3^1 + e9 



27r2 



120 



In the same manner, the expression for the pressure can be found to be 



30 



7 Ktt 
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{Pu + Pd) 



(9) 



(10) 



In the non-relativistic limit, on the other hand, we shall expand e as e ~ am(l 



Unfortunately, we do not have any expression for the even moments of degenerate distribu- 
tion function in terms of a finite power series of ^. Therefore we keep the integral expressions 
and write the energy density as 



Pv + Pv 

and the pressure as 
3(P. + Pp) 



27r2 



am j q^dq{f^ + /p) + j q'^dq{f^ + /p) 



{pvo + Ppo) iamQ2 + 77^Q4 ) , 
\ 2am J 



(11) 



{ksTu 



27r2 



— / q'dqiU + U) - 
am J 

{Pvo + Ppo) ( — ©4 - 777— T^0f 

\am 2[amy 



2{amy 



q'dq{U + U) 



(12) 
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FIG. 1: The neutrino mass as a function of ^ normalized by ^ = non-degenerate neutrino. With 
^^h'^ fixed, the mass of neutrinos in lepton asymmetric models are lighter than that in symmetric 
models. 



Here we have defined the even-moments of distribution function 6„, normalized by the third 



moment of non-degenerate fermion distribution, J q^dq-^ ^ 



(l+e9) 



120 



(13) 



Before finishing this subsection we would like to comment on the relation between the 
cosmological neutrino density parameter Vl^ and the mass of a neutrino mj, 2l|. Basically, 
as long as m^, is larger than ksTy, the neutrino mass is related to VL^, as 



rrty 



(14) 



where VLy is the density parameter of neutrinos, pcrit is the critical density, Ny is the number 
of species of massive neutrino, ny is the number density of one species of massive neutrino. In 
the case of degenerate neutrino, the number density is different from that of non-degenerate 
one, and it is given by 



27r2 



X 2 



d^q 
(2^ 



iUq) + IM) 



120 



(15) 



By remembering that for non-degenerate massive neutrinos the number density is expressed 

as 

, , , {kBTuf „ 3 



, X 2 X -C(3) 
27r2 2^^ ' 



(16) 



where C, is the zeta function, we obtain the following relation between masses of degenerate 
neutrino and non-degenerate one as 



m,. 



|C(3) 
1^02(0 



(17) 
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Here ^ m^g = rni,^^ + m^^^+m^^^^ = 94.1 eV Q^h'^ is the neutrino mass in the non-degenerate 
cases, and Ty(0) and T,^{^) are the present temperatures of standard and degenerate neutri- 
nos, respectively. 

The temperatures of degenerate and non-degenerate neutrinos can differ greatly from 
each other if the chemical potential of neutrinos is so large that neutrinos decouple from the 
rest of the plasma before muon-antimuon annihilation. In reality, however, neutrinos can 
be kept in kinetic equilibrium with even with a significant degeneracy to hold the same 
temperature as the rest of the plasma until at least the muon-antimuon annihilation ends 



|22l |. while number freeze-out can occur before cosmological QCD phase transition [23|, |24 

1 /3 

So, we can safely set T,^(0) = T,^(0 = (^) ^7,0 in Eq.(|T71). Therefore, with the density 
parameter Qi^h'^ being fixed, the mass should be monotonically smaller as the neutrinos are 
degenerate larger. The mass ratio between degenerate and non-degenerate cases is plotted 
in Fig. [H 



B. perturbations 

In this paper we only consider the tensor type perturbations. The line element is given 

by 

ds^ = gf.^dx^dx" = a^(r) [-dr'^ + {6ij + hjj)dx'dx^] , (18) 

where hjj denotes tensor type perturbation (gravitational waves) around flat Friedmann- 
Robert son- Walker metric. We shall work in the transverse traceless gauge, which is set by 
the conditions /i^- j = hj^ = 0. The evolution of gravitational waves hj^ is determined by the 
linearized Einstein equations: 

hi + 2 Q hi + k^hl = WTcGa^TTij , (19) 

where dots indicate conformal time derivatives, and Hij is the anisotropic part of the stress 
tensor Tij, defined by 

a^TTij = Tij - pQij . (20) 

Following the standard procedure, we expand the distribution function of neutrinos into 
homogeneous and perturbed inhomogeneous parts as 

fu + f, = {fu + ff>)il + ^^ix\T,q,)) . (21) 

Note that we take comoving momentum in the tetrad frame qi as our momentum variable, 
which is related to the (spatial part of) four momentum as qi = a?P\5ij + \hij). By 
virtue of this momentum variable, one can express the moments of distribution function such 
as energy density, pressure and anisotropic stress in a simple manner, i.e., qi can be treated 
as that on the flat space background. For example, the anisotropic stress of neutrinos can 
be written as 

= {kT.r I 7(7.7, - + (22) 

where we write the comoving momentum qi in terms of its magnitude and direction: qi = qji, 
with 5*-'7i7j = 1. 
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The perturbed distribution function satisfies tensor type collisionless Boltzmann equa- 
tion in fc-space 

where is a directional cosine between k and q, and logarithmic derivative of the distribution 
function with respect to q is given by 



d ln(/j, + fp) q A 1 + cosh ^ cosh q 



dliaq e~'i + cosh ^ \ cosh ^ + cosh q 

One finds the solution of Eq. (123!) as 



(24) 



2 rflng 

where new time variable u is introduced by 



= + / rf«'e-(-'-)/.:,(«0 , (25) 



u = u{q) = k [ . (26) 



One of the effects of finite mass comes through Eq.( !26l) . which delays the timing of new time 
variable u compared with the conformal time kr. 

By inserting the solution Eq. (!25l) into Eq. (!22l) . and integrating by parts, one obtains 



TT, 



(27) 

One can easily confirm that, in cases of massless non-degenerate neutrinos, Eq. (l271) reduces 
to the results obtained in [15|, ll7|. By substituting this expression into the l.h.s. of Eq. (|T9il . 
one sees that the evolution of gravitational waves is given by solving an integro-differential 
equation. 

In order to see how neutrino asymmetry affects the evolution of gravitational waves, let 
us first consider massless degenerate neutrinos. When m = 0, the expression for nfj reduces 
to 

■5,n.assless = '^P^ + P.) £ du' ^^^^^^^^^^""'''^ ' (28) 

In such cases, therefore, q dependence of the distribution functions (Eq. (^^) can be 
integrated away into the energy density, just like in the non-degenerate cases. So the effect 
of ^ arises only through the background energy density (in the first term of Eq. (jHD), 
and is completely described by introducing an effective number of massless non-degenerate 
neutrinos. 

For massive degenerate neutrinos, however, the situation is quite different. One can no 
longer express the momentum integration in Eq. fl27|) in terms of well-known moments such 
as density, pressure and so on. In other words, it will not be possible to renormalize the 
dependence of ^ by the effective number of massless non-degenerate neutrinos. Thus, in 
principle, these two contributions are distinguishable. 
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In the relativistic limit, the expression of nf, can be expanded to 



T 

TT, ^ 



a?rn? 



TT 



X / dn' / ' c^xe"("-"') f^^ l.h,,{u') . (29) 



16 9ti 

In the non-relativistic limit, on the other hand, vr^^- can be expressed as 



J -I 

T 



For both limits, a finite mass always acts to suppress the anisotropic stress relative to the 



massless case. One finds that when neutrinos become massive, g/e — > 0, then tiJ- 0. 



III. EVOLUTION OF GRAVITATIONAL WAVES 

In the absence of neutrino free streaming in which vr^ = 0, the evolution equation of 
gravitational waves (Eq. (fT9|) ) with an appropriate initial condition has well-known analyti- 
cal solutions: h oc joikr) in the radiation dominated era, and h oc ji{kT) / [kr) in the matter 
dominated era. Thus for both eras, the amplitude of gravitational waves diminishes as 
with oscillations when waves are well inside the cosmic horizon [kr ^ 1). This is a general 
feature of waves in an expanding universe. 

The main effect of neutrinos on this evolution of gravitational waves is the damping 
of amplitude by their free streaming |[17|], and it can be significant when the universe is 
radiation dominated and neutrinos are relativistic. Because the gravitational waves with 
wavenumber k > 0.01 Mpc~^ enter the Hubble horizon before matter-radiation equality 
for the case of standard massless neutrinos, they suffer from significant damping. This 
condition holds true unless mass of neutrinos is heavier than ~ 1 eV, in which cases neutrinos 
become non-relativistic before the equality. On the contrary, the gravitational waves with 
longer wavenumber k < 0.01 Mpc~^ are not significantly affected by neutrinos. Examples 
of evolutions of gravitational waves with wavenumbers k = 0.005, 0.5 Mpc~^ are illustrated 
in Fig. H 

In Fig. [21 we also show the effect of non-zero chemical potential on the evolution of 
gravitational waves. Non-zero chemical potential modifies the evolution of gravitational 
waves mainly in two ways through an effective increase of energy density of neutrinos. One 
is an effective increase of anisotropic stress in Eq. (128|) . This increases the effect of free 
streaming neutrinos, and consequently, the damping of gravitational waves becomes larger. 
The effect is clearly seen in the numerical evolution of gravitational waves (black dash- 
dotted line) in Fig. [2J The other is the shift of matter-radiation equality to the later time 
in the history of the expanding universe. This indirectly gives rise to an overall increase in 
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FIG. 2: The effects of neutrino free streaming and ^ in the time evolution of gravitational waves. 
The black lines show the numerical solutions of gravitational waves with wavenumber A; = 0.5 
Mpc~^, and blue (gray) lines show those with wavenumber k = 0.005 Mpc~^. The two thickest 
lines are solutions without neutrinos (no v). The dash-dotted lines are numerical solutions with 
degeneracy parameter ^ = 2. The straight thick and dashed lines are those without chemical 
potentials. The lines with smaller wavenumber (k = 0.005) are almost indistinguishable in this 
figure. 



amplitude of gravitational waves at present for modes which have come across the horizon 
in the radiation dominated universe. 

Let us quantitatively discuss the former effect first. An analytic estimation of the damping 
of gravitational waves by neutrino free streaming was given by [1^ in a sophisticated way. 
Here we follow their argument and apply it to the lepton asymmetric cosmology. They 
expanded a solution of gravitational waves h{u) in a series of spherical Bessel functions as 

oo 

h{u) = ^anjn{u) . (31) 

n=0 

At sufficiently late times, the amplitude of gravitational waves asymptotically takes of the 
form 

h{u) = Asm{u + 6)/u (32) 

For large argument, all of the even order Bessel functions go as ±smx/x so the damping 
factor A is given by 

5 

A ^ 5^(-l)'^a2. . (33) 

n=0 

Coefficients a2n are determined by the recursion relation 

1 ar J2k=0 0'2kC2k,2n-2k 
a2n - -i.O/,,-— — - -. , (^4j 

2n{2n + 1) + l-6ju 

starting with oq = 1, where f^, is the fraction of the energy density in neutrinos. Numerical 
coefficients C2k,2n-2k are listed in [2^. 

For the standard (massless and non-degenerate) neutrinos, = 0.40523, which gives 
A = 0.80313 251]. However, in the lepton asymmetric cosmology, it is not true in general. 
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FIG. 3: The damping factor as a function of degeneracy parameter The damping effect by 
neutrino free streaming becomes larger with larger ^. 

Using degeneracy parameter ^, the fraction of the energy density in neutrinos is now given 
by 

^•"^^^ ^ F(0 + 1.46773 ' ^^^^ 

where 

When ^ = 7r, we get A ^ 0.64231. This enhancement of damping by neutrino free streaming 
is illustrated in Fig. [31 

Next, we discuss the second effect, i.e., the shift of matter-radiation equahty. This can 
result in a modification of the spectrum of gravitational waves at present. In order to 
investigate it, let us introduce the power spectrum of gravitational waves as 

Al^^dM"). (37) 

In the simplest inflationary models, it is expected that the power spectrum is nearly scale 
invariant and has the amplitude of ^ (^^^1^ ' ^'^^re i/inf is the Hubble parameter during 

inflation. Alternatively, the spectral energy density, flh{,k) oc k^A\, is often used instead of 
the power spectrum. It has been shown that, as a result of the evolution of gravitational 
waves from the deep radiation dominated era through the matter dominated era to the 
present universe, the spectral energy density has a fiat spectrum at > 0.01 Mpc~^, where 
the corresponding Fourier modes have come across the cosmic horizon during radiation 
dominated era, and at k < 0.01 Mpc^^, where those have crossed the horizon during 
matter dominated era [1^. We shall discuss the spectrum of gravitational waves in terms 
of this spectral energy density. 

Because the energy density spectrum of gravitational waves, whose corresponding Fourier 
modes cross the cosmic horizon during the matter dominated epoch, has a slope of if the 
primordial spectrum was scale invariant, the amplitude will be boosted by a shift of matter- 
radiation equality, provided that the normalization of overall amplitude of gravitational 
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FIG. 4: Spectra of gravitational waves with and without lepton asymmetry. The net effect from 
lepton asymmetry is a rise of the plateau at A; > 10^^ Mpc^^. 



waves is given at superhorizon scales. The amplification factor is therefore given by, 

_[ ^ _ 2 + 1.36F(0 



^h, standard \ ^eq,standard / 3.36 



(38) 



To derive this relation we have used the fact that the redshift at the matter-radiation equality 
is given by 

^ , _ Pcrit^m _ 8.83 X 10=^ 

P.,0 + P..0 2 + 1.36F(0 V 0.11 ) ' ^^^^ 
and the corresponding comoving wavenumber is therefore 

= = 0.0145 [Mpc-i] ( ^^^^1 ( ^1 . (40) 

^ 1 + Zeq ^V2 + l-36F(Oy V 0.11 / ^ 

Here flm is total matter density in units of the critical density pcrit, and p^^o, Pfo are en- 
ergy densities of photons and neutrinos today. When ^ = vr, ^ 3.6. This indirect 

amplification of gravitational waves dominates the damping effect discussed earlier. This is 
illustrated in Fig. HI 

As mentioned earlier, if the mass of neutrinos are heavier than ~ 1 eV, neutrinos have 
gone out from the ultra-relativistic regime before the matter radiation equality and then 
made their free streaming effect on gravitational waves smaller. In Fig. [5l we depict the 
effect of finite mass of neutrinos on the time evolution of gravitational waves. From a 
practical viewpoint, however, the mass of neutrinos as large as > 1 eV is not realistic because 
even CMB anisotropy data alone by WMAP have already put constraint on the mass as 
< 2 eV (at 95% confidence) [23] and the constraint can be tighter ^rrii, < 0.68 
eV if the information from the matter power spectra is included Therefore we may 
conclude that the effect of finite mass of neutrinos cannot be significant on the evolution of 
gravitational waves. 



11 




FIG. 5: The effects of finite neutrino mass {nii, = 4 eV) in the evolution of gravitational waves with 
the wavenumber k = 0.05 Mpc^^. The finite mass suppresses the anisotropic stress of neutrinos 
and the damping effect becomes smaller. 



IV. DISCUSSION AND SUMMARY 

In this paper we investigated the evolution of gravitational wave background in lepton 
asymmetric cosmology. We showed that the non-zero degeneracy of neutrinos can modify 
the present spectrum of gravitational wave background through two ways. One is that the 
non- negligible chemical potential results in the larger anisotropic stress of neutrinos and 
thus larger damping of amplitude of gravitational waves. The other is the shift of matter 
radiation equality, which indirectly raises the plateau of the spectrum at A; > 10^^ Mpc^^. 
These effects can act in the opposite direction from each other. In fact, we showed that the 
latter effect dominates the former, and the spectrum has a larger power in total than that 
in the case with non-degenerate neutrinos. These effects follow from the effective increase 
of number of massless neutrinos. We presented straightforward numerical calculations and 
analytic estimates for both of the effects. 

An important implication of our results is that, should the other cosmological parameters 
such as Vtmh'^i Ut, m,^ ever be known, one can determine in principle the lepton asymmetry 
by measuring the gravitational wave amplitude around k > 10~^Mpc~^. This fact would 
become more hopeful noting that, among these cosmological parameters, neutrino mass 
would not be so important because the mass itself does not affect the shape of the spectrum 
of gravitational waves once it is constrained to be less than ~ eV from the other cosmological 
or terrestrial experiments. 

Of course, the detection of gravitational wave background is promising but still very 
challenging. Because the damping due to neutrino free streaming can be found only at 



lower frequency than ~ 2 x 10~^° Hz 15| assuming standard decoupling of neutrinos, it 
will be difficult to detect it in a direct way. Indeed, although the critical frequency would 
be higher as a consequence of smaller Hubble horizon scale at neutrino decoupling in the 
presence of a significant degeneracy, the frequency will be still far lower than that aimed by 
the proposed gravitational wave experiments. The most promising way to measure the lepton 
asymmetry would be to detect the gravitational wave background indirectly by observing 
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the curl modes of cosmic microwave background polarization patterns. It is expected that, 
in the future planned experiments of CMB anisotropics such as Planck or CMBpol, the 
lepton asymmetry would be determined if the degeneracy parameter is as large as ^ ~ 0.1 
and the contamination from divergence mode of polarization by cosmic shear is cleaned out. 

As stated in the introduction, in case that flavor equilibrium between all active neutrino 
species is established before BBN epoch, which is indicated by the large mixing angle solution 
of the solar neutrino problem, light element abundances have already put a tight conservative 
limit on degeneracy parameter as |^| < 0.1 [ll[. However, because the other solution may 
exist in which flavor mixing can be achieved only partially [ll|, and because the evolutions of 
CMB anisotropics are completely independent from BBN physics, it will still be important to 
investigate how much the sensitivity of CMB anisotropics for ^ can be improved by including 
the effect of on the evolution of gravitational waves studied here. 

Acknowledgments 

JY thank Galileo Galilei Institute for Theoretical Physics for the hospitality and the 
INFN for partial support during the completion of this work. This work was partially 
supported by JSPS Grant-in-Aid for Scientific Research Nos. 1809940(KI), 18740157(MY), 
and 16340076 (JY). M.Y. is supported in part by the project of the Research Institute of 
Aoyama Gakuin University. 



[1] For review, see E. W. Kolb and M. S. Turner, The Early Universe (Addison- Wesley, New 
York, 1990); K. A. Olive, G. Steigman, and T. P. Walker, Phys. Rep. 333, 389 (2000); D. 
Tytler, J. M. O'Meara, N. Suzuki, and D. Lubin, Phys. Scr., T85, 12 (2000). 

[2] C. L. Bennett et al, Astrophys. J., Suppl. Ser. 148, 1 (2003);148, 175 (2003); H. V. Peiris et 
al, Astrophys. J., Suppl. Ser. 148, 213 (2003). 

[3] D. Spergel et al, astro-ph/0603449. 

[4] S. Yu. Khlebnikov and M. F. Shaposhnikov, Nucl. Phys. B308, 885 (1988); J. A. Harvey and 

M. S. Turner, Phys. Rev. D 42, 3344 (1990). 
[5] K. Enqvist, K. Kainulainen, and J. Maalampi, Phys. Lett. B 244, 186 (1990); R. Foot, M. J. 

Thomson, and R. R. Volkas, Phys. Rev. D 53, 5349 (1996); X. Shi, ibid. 54, 2753 (1996); For 

review and further references, see A. D. Dolgov, Phys. Rep. 370, 333 (2002), Sec. 12.5. 
[6] J. A. Harvey and E. W. Kolb, Phys. Rev. D 24, 2090 (1981); A. Casas, W. Y. Cheng, and G. 

Gelmini, Nucl. Phys. B538, 297 (1999); J. McDonald, Phys. Rev. Lett. 84, 4798 (2000). 
[7] J. March-Russell, A. Riotto, and H. Murayama, J. High Energy Phys. 11, 015 (1999). 
[8] M. Kawasaki, F. Takahashi, and M. Yamaguchi, Phys. Rev. D 66, 043516 (2002). 
[9] M. Yamaguchi, Phys. Rev. D 68, 063507 (2003); T. Chiba, F. Takahashi, and M. Yamaguchi, 

Phys. Rev. Lett. 92, 011301 (2004); F. Takahashi and M. Yamaguchi, Phys. Rev. D 69, 083506 

(2004). 

[10] S.H. Hansen, G. Mangano, A. Melchiorri, G. Miele, and O. Pisanti, Phys. Rev. D 65, 023511 
(2002). 

[11] A.D. Dolgov, S.H. Hansen, S. Pastor, S.T. Petcov, G.G. Raffelt, D.V. Semikoz, Nucl. Phys. 
B632, 363 (2002). 



13 



[12] K. S. Babu, Phys. Lett. B 275, 112 (1992); L. Bento and Z. Berezhiani, Phys. Rev. D 64, 

115015 (2001); A. D. Dolgov and F. Takahashi, Nucl. Phys. B688, 189 (2004); 
[13] A. A. Starobinsky, JETP Lett. 30, 682 (1979). 
[14] N. Seto and J. Yokoyama, J. Phys. Soc. Japan, 72, 3082 (2003). 
[15] Y. Watanabe and E. Komatsu, Phys. Rev. D 73, 123515 (2006). 
[16] D. J. Schwarz, Mod. Phys. Lett. A 13, 2771 (1998); Ann. Phys. 12, 220 (2003) 
[17] S. Weinberg, Phys. Rev. D 69, 023503 (2004). 
[18] J. Lesgourgues and S. Pastor, Phys. Rev. D 60, 103521 (1999). 

[19] M. Orito, T. Kajino, G. J. Mathews and Y. Wang, Phys. Rev. D 65, 123504 (2002). 

[20] M. Lattanzi, R. Ruffini, and G. V. Vereshchagin, Phys. Rev. D 72, 063003 (2005). 

[21] K. Preese, E. W. Kolb, and M. S. Turner, Phys. Rev. D 27, 1689 (1983). 

[22] K. Ichikawa and M. Kawasaki, Phys. Rev. D 67, 063510 (2003) 

[23] H.-S. Kang and G. Steigman, Nucl. Phys. B372, 494 (1992) 

[24] M. Orito, T. Kajino, G. J. Mathews and R. N. Boyd, astro- ph/0005446^ 

[25] D. A. Dicus and W. W. Repko, Phys. Rev. D 72, 088302 (2005). 

[26] B. Allen, Some Topics on General Relativity and Gravitational Radiation. Proceedings of the 
"Spanish Relativity Meeting '96", Valencia, Spain, September 10-13, 1996, edited by J. A. 
Miralles, J. A. Morales, and D. Saez, 1997., p. 3. 

[27] M. Fukugita, K. Ichikawa, M. Kawasaki, O. Lahav, Phys. Rev. D 74, 027302 (2006). 



14 



